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The global burden of stroke continues to grow. Although stroke prevention strategies (e.g., 
medications, diet, and exercise) can contribute to risk reduction, options for acute interven- 
tions (e.g., thrombolytic therapy for ischemic stroke) are limited to the minority of patients. 
The remaining patients are often left with profound neurological disabilities that substan- 
tially impact quality of life, economic productivity, and increase caregiver burden. In the last 
decade, however, the future outlook for such patients has been tempered by movement 
toward the view that the brain is capable of reorganizing after injury. Many now view brain 
recovery after stroke as an area of scientific research with large potential for therapeutic 
advances, far into the future (Broderick and William, 2004). As a probe of brain anatomy, 
function and physiology, magnetic resonance imaging (MRI) is a non-invasive and highly 
versatile modality that promises to play a particularly important role in such research. Here 
we provide a basic review of MRI physical principles and applications for assessing stroke, 
looking toward the future role MRI may play in improving stroke rehabilitation methods and 
stroke recovery. 
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INTRODUCTION 

Diagnostic imaging is an invaluable aspect of clinical stroke med- 
icine, providing the location, volume, and the nature of the stroke 
lesion. Anatomical imaging of stroke patients historically has been 
done with computed tomography (CT), but ever-increasingly is 
supplanted by the superb soft tissue contrast provided by magnetic 
resonance imaging (MRI). The versatility of MRI methods also 
enables much more detailed biophysical information on stroke 
physiology, above and beyond lesion structure. For example, just 
after stroke onset, diffusion-weighted MRI (DWI), and perfusion- 
weighted MRI (PWI) provide information about diffusion of 
water molecules and microvascular blood flow within brain tis- 
sue, respectively. The DWI and PWI methods help to evaluate 
the ischemic zone surrounding infarcted tissue that is poten- 
tially salvageable by recanalization therapies. Magnetic resonance 
angiography (MRA) approaches are also available to character- 
ize larger-scale vasculature. In the post-acute and chronic phases, 
functional MRI (fMRI) offers the ability to detect alterations in 
brain activation patterns post stroke, either associated with a par- 
ticular behavioral task or during the resting state. It is also possible 
now to image cerebral blood flow completely non-invasively using 
arterial spin labeling (ASL) MRI. In this review, we provide bio- 
physical understanding of these and other basic MRI methods, 
and discuss their application to stroke recovery. 

BASIC PHYSICAL PRINCIPLES 

Use of MRI in medicine has evolved from the discovery of nuclear 
magnetic resonance (NMR) in the 1950s (Bloch, 1946; Purcell 
et al., 1946). The term NMR is apt- nuclear refers to specific atomic 



nuclei which, when placed in a large applied static magnetic field, 
can undergo resonance (substantial energy absorption and subse- 
quent energy emission if a pulsed magnetic wave is applied at a 
certain characteristic frequency). For NMR to occur, the nuclei of 
interest must have the quantum mechanical property of non-zero 
"spin." Such nuclei simplistically can be thought of as very small 
objects with spinning charge, and therefore exhibiting nuclear 
magnetic fields. When placed in a large applied magnetic field, 
the nuclear magnetic fields tend to align with the applied field. 
The simplest scenario involves the hydrogen nucleus, or proton, 
which can take on two different energy states: a low-energy "spin- 
up" state in which the nuclear magnetic field partially aligns to the 
external field, and a higher-energy "spin-down" state in which the 
nuclear field aligns partially anti-parallel. The word "partial" refers 
to the fact that in each state, the alignment is not perfect. Rather, 
the nuclear field "precesses" around the direction of external field. 
A good analogy for this effect is a spinning top, the axis of which 
begins to wobble (precess) about the direction of gravity as the top 
"winds down." In the case of NMR, the frequency of precession 
is known as the Larmor frequency. The Larmor frequency is the 
resonance frequency of the system, and transitions between energy 
states are possible if a magnetic pulse is applied at the Larmor fre- 
quency. Typically, a pulse in the radiofrequency (RF) portion of 
the electromagnetic spectrum is required for NMR, similar to the 
frequencies used in FM radio. 

The NMR effect is one remarkable example of how quan- 
tum mechanical effects occurring over extremely small distances 
(approximately 1CP 10 m) can manifest as macroscopic, observable 
signals. Within a material containing NMR-capable nuclei, the 
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nuclear magnetic fields for the given energy states can be thought 
of as "adding up," resulting in a net nuclear magnetization. The 
material is very slightly magnetized when placed in a large static 
magnetic field, due to the population imbalance in favor of pro- 
tons in the more favorable low-energy state. The magnetization is 
typically small at room temperature because thermal energy is suf- 
ficient to promote a small change between the spin energy states. 
If the external magnetic field is made strong enough, the nuclear 
magnetization increases to the point that a signal can be easily 
observed above the noise background of a detector circuit. For 
this reason, current MRI systems operate using superconductor 
technology to generate very large static magnetic fields of 1.5-3.0 
Tesla, with ultra-high field systems at 7.0 Tesla and beyond under 
development (1 Tesla is approximately 2 x 10 4 times the Earth's 
magnetic field). 

The conventional manner for observing nuclear magnetization 
involves three steps. First, the sample is placed in a large applied 
magnetic field, creating magnetization in the parallel or "longi- 
tudinal" direction. Second, a magnetic RF pulse is applied at the 
Larmor frequency, causing energy absorption that perturbs the 
energy states and consequently the magnetization. The RF pulse 
has the macroscopic effect of "tipping" the magnetization out of 
alignment with the applied field, leading to a vector component 
oriented transverse (orthogonal) to the longitudinal direction and 
that oscillates at the Larmor frequency. Lastly, the oscillating trans- 
verse magnetization component is detected as an NMR signal 
using a "receiver coil." In its simplest form, a receiver coil consists 
of an appropriately oriented loop of wire that records an NMR- 
related oscillating voltage induced by the transverse magnetization, 
according to Faraday's Law of Induction. 

Early in the development of NMR, it was recognized that hydro- 
gen exhibits the largest nuclear magnetic field and thus the largest 
NMR signal of all the elements. This is fortuitous, given that hydro- 
gen in water is by far the most biologically abundant of nuclei. An 
additional strength is the remarkable versatility of NMR, to probe 
not only the magnetic field properties of the nucleus, but also mol- 
ecular interactions. Collectively, these effects provide an unrivaled 
mechanism to probe the soft tissues of the body in a non-invasive 
manner. 

MAIN SOURCES OF NMR SIGNAL CONTRAST 

In the early 1970s, it became increasingly apparent that biological 
tissues often exhibit markedly different NMR signal characteris- 
tics (Damadian, 1971). There are several ways to obtain tissue 
contrast, i.e., the ability to detect different tissue types based on 
the associated NMR signal differences. First, the number of pro- 
tons per unit volume, or proton density, is a factor that weights 
all NMR signals. Typically the proton density of tissues is rather 
uniform throughout the body and a moderate source of signal 
contrast. Much improved signal contrast is obtained based on sev- 
eral NMR parameters that characterize how magnetization tends 
to recover toward its equilibrium condition after the RF pulse is 
turned off. The RF pulse tips magnetization out of alignment with 
the external magnetic field, typically producing vector compo- 
nents in the longitudinal direction and in the transverse direction. 
Magnetization then immediately starts to return to its equilibrium 
condition, resulting in recovery of longitudinal magnetization, and 



decay of transverse magnetization. In simple fluids such as water, 
both processes are exponential in nature. Recovery of longitudinal 
magnetization is characterized by the time constant Tl, whereas 
decay of transverse magnetization is characterized by the time con- 
stant T2. In tissues, typical Tl values range from approximately 
300-2000 ms at 1.5 and 3.0 Tesla, respectively, with a tendency 
for Tl to increase as the applied magnetic field increases. The 
analogous values for T2 range from approximately 30-300 ms. 
Table 1 provides proton density, Tl, and T2 values for different 
brain tissues. 

Table 1 indicates that Tl and T2 values differ by approximately 
a factor of 10, and vary for different tissues. The reason for tissue- 
related differences in relaxation times is multifaceted (Bronskill 
and Graham, 1993), but several of the key issues can be briefly 
summarized. Considering Tl relaxation first, the recovery of lon- 
gitudinal magnetization is primarily influenced by re-population 
of spin states. When an RF pulse causes increased population of the 
spin-down state, the subsequent probability of spontaneous decays 
to the low-energy spin-up state is very low. Transitions occur pre- 
dominantly due to "stimulated emission," an effect that requires 
a source of magnetic field fluctuations primarily at the Larmor 
frequency. The source of these fluctuations is dynamic proton- 
proton magnetic field interactions from within water molecules 
and between neighboring water molecules as they tumble and 
translate. This relaxation process is influenced by factors such as 
temperature, viscosity, and the presence of ions, macromolecules, 
and tissue microstructures such as cell membranes. The effect is 
additionally complicated by the fact that water molecules typically 
can diffuse through multiple different relaxation environments 
within tissues during the Tl timescale. The process of T2 relax- 
ation is additionally influenced by the magnetic fields associated 
with very slow components of molecular motion, as these perturb 
the Larmor frequency for individual protons and reduce the coher- 
ence of nuclear magnetic fields in the transverse plane. This effect 
explains why T2 is much shorter than Tl in tissues, and why pro- 
tons from large macromolecules are typically not observed during 
MRI. In the latter case, these protons have very slow rates of mol- 
ecular motion and the resulting magnetic field non-uniformity at 
the nuclear scale leads to T2 values of approximately 10 (is, too 
small for observation with conventional MRI technology. 



Table 1 | List of relaxation times by tissue type and main magnetic 
field svtrength (Gati et al., 1997; Kruger et al., 2001; Lu et al., 2005; Wu 
and Wong, 2006). 

Field Tissue T1 (ms) T2 (ms) T2* (ms) Proton 

strength (T) density 



White matter 


510 


67 


78 


0.61 


Gray matter 


760 


77 


69 


0.69 


Arterial blood 


1441 


290 


55 


0.72 


CSF 


2650 


280 


n.a. 


1.0 


White matter 


1080 


70 


50 


0.61 


Gray matter 


1820 


100 


50 


0.69 


Arterial blood 


1932 


275 


46 


0.72 


CSF 


3817 


1442 


n.a. 


1.0 
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Table 1 also includes one other NMR parameter, known as 
T2*. Both T2 and T2* parameters reflect decay of transverse 
magnetization, as measured using different NMR experiments. 
T2 is measured using a technique known as "spin-echo refo- 
cusing" (Carr and Purcell, 1954) that measures the decay of 
transverse magnetization while compensating for all macroscopic 
sources of magnetic field inhomogeneity in space (e.g., spatial 
variation in the applied magnetic field). The T2* parameter is 
measured in the absence of such compensation, and is typi- 
cally smaller than the T2 value. The relatively fast T2* decay 
is due primarily to the fact that biological tissues are highly 
heterogeneous with the constitutive components often support- 
ing slightly different applied magnetic fields at the microscopic 
level. These differences arise from spatial variations in mag- 
netic susceptibility, a property of all materials. The result is 
that tissue will exhibit a range of different Larmor frequencies 
even if the external main magnetic field is perfectly uniform, 
reducing the coherence of signal oscillation in the transverse 
plane and causing more rapid decay of transverse magnetiza- 
tion than is measured by the T2 parameter. Thus, T2* provides 
an additional mechanism to achieve tissue-specific NMR signal 
contrast. 

DEVELOPMENT OF MRI 

Early investigations of NMR signals from tissues focused partic- 
ularly on the ability to detect cancer and distinguish malignant 
from benign lesions, and provided some of the impetus for devel- 
oping methods to encode NMR signals spatially in the form of 
images. Key innovators in these subsequent developments were 
Lauterbur and Mansfield, who received the Nobel Prize in phys- 
iology and medicine in 2003. Over the years, the term "nuclear" 
has been dropped to avoid any potential misconception that the 
imaging methodology involves use of ionizing radiation - hence 
the modern acronym MRI. The first whole-body MRI systems 
were developed for commercial use in the late 1980s, and since 
then there has been a steady increase in technical capabilities. 
Today, multifaceted MRI "protocols" are achieved with approx- 
imately 1 mm spatial resolution and exquisite image contrast. 
Specialized sequences of RF pulses are used to provide images 
weighted by different MRI signal parameters in examination times 
of approximately 30 min. 

The key technological concept that led to the development of 
MRI is the "imaging gradient" for spatial encoding. The imag- 
ing gradient is an electrical circuit that produces a linear change in 
magnetic field as a function of position. MRI systems provide mag- 
netic field gradients in the three orthogonal directions x, y, and z. 
Importantly, although the gradient directions are orthogonal, the 
orientation of the magnetic field from each gradient is always in 
the direction of the large external magnetic field. Consequently, 
a given gradient produces a linear variation in Larmor frequency 
with position and can be used to encode NMR signals in the gradi- 
ent direction (Lauterbur, 1973). In comparison with the external 
magnetic field, gradient fields are much smaller, on the order of 
lOmilliTesla/meter. They are also designed to be pulsed on and 
off throughout the imaging procedure, enabling spatial encoding 
in the x, y, and z directions. It is the pulsing of imaging gradi- 
ents (leading to vibrational forces within the magnet bore) that 



is responsible for the characteristic buzzing, knocking, or pinging 
sounds that occur during MRI scans. 

Both "multi-slice" and three-dimensional (3D) MRI require 
multiple repetitions of RF pulse excitation, data acquisition, and 
recovery periods to obtain all the necessary spatial encoding infor- 
mation required for image reconstruction. Typically, 3D MRI is 
advantageous over 2D MRI because isotropic spatial resolution is 
possible and the data can be straightforwardly reformatted post- 
acquisition into different image plane orientations (e.g., axial, 
sagittal, coronal, arbitrary oblique). However, 2D MRI is typically 
more time- efficient and is widely used on this basis. Clinically, 
both types of imaging have scan times that usually fall within the 1- 
10 min range. There is also a class of "fast" MRI techniques capable 
of providing a fully encoded image slice in approximately 100 ms, 
or multi-slice datasets in approximately 1-2 s. The most common 
of these techniques is echo planar imaging (EPI) (Stehling et al, 
1991), capable of full spatial encoding within a single RF excita- 
tion followed by data acquisition during rapidly varying gradient 
waveforms. EPI can be used to image time-dependent physiolog- 
ical processes that occur on the timescale of seconds, such as the 
wash-in of injectable contrast agents. However, fast imaging has 
costs: reduced signal-to-noise ratio (SNR), reduced spatial resolu- 
tion, and regions of signal loss and substantial geometric distortion 
arising from sensitivity to magnetic field uniformity. 

MRI OF CEREBROVASCULAR RISK FACTORS 

Magnetic resonance imaging is being used increasingly as part of 
the cerebrovascular work-up to assess risk of stroke. Vascular risk 
factors like hypertension, type 2 diabetes mellitus, and the natural 
course of aging are associated with an increase in the prevalence of 
subclinical brain lesions (Thompson and Hakim, 2009). For exam- 
ple, individuals with type 2 diabetes that have brain lesions on MRI 
appear to have a two to fivefold increase in stroke occurrence com- 
pared to non-diabetics (Baird et al, 2002). Cerebrovascular risk 
factors can manifest as white matter disease (see Figures 1E,F), 
lacunes, small vessel disease, and/or Virchov Robin spaces, which 
may be considered subclinical if they do not produce an overt neu- 
rological symptom. These features may be viewed as an "incidental 
finding" but it is more apt to refer to these events as "covert strokes." 
For example, meta-analyses show that white matter lesions are a 
significant risk factor for stroke, dementia, and death (Biernaskie 
et al., 2004). Thus, as we advance our understanding of the stroke 
in the acute and chronic stages it is imperative to characterize 
the underlying cerebrovascular pathologies that may have con- 
tributed to the stroke in the first place, and that will also influence 
the potential for recovery. 

MRI OF ACUTE STROKE 

An MRI protocol is a set of individual image acquisitions that 
each have different signal contrast weightings, assembled with the 
goal of detecting pathology in a comprehensive manner. Such 
protocols are becoming increasingly critical in an acute stroke 
setting, to direct stroke interventions and characterize recov- 
ery (Leiva-Salinas et al., 2011). For example, it is possible to 
design an MRI protocol that identifies regions of possible active 
ischemia, hemorrhage, and occluded vessels, and that provides 
insight regarding the size and location of the core infarct as well as 
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FIGURE 1 | Multiple MRI pulse sequences can be used to assess 
cerebrovascular risks. The magnitude (A) and phase (B) images in phase 
contrast angiography can be used to calculate cerebral blood flow velocity in 
arteries (red arrows) and veins. (C) Susceptibility-weighted images (SWI) 
provide high spatial resolution (<1 mm in plane resolution) of venous 



the extent of surrounding tissue that potentially can be salvaged 
(the ischemic penumbra). A typical acute stroke MRI protocol 
consists of diffusion-weighted imaging (DWI), fluid attenuated 
inversion recovery (FLAIR), MRA, perfusion-weighted imaging 
(PWI), and T2-weighted imaging. Some representative examples 
of these techniques are shown in Figures 1 and 2, as they are sen- 
sitive both to overt and covert stroke, with further detail to follow 
below. 

Standard MRI sequences that provide Tl-weighting or T2- 
weighting, and their derivatives, are well known to be insensitive 
to the immediate effects of cerebral ischemia. By contrast, DWI 
is highly sensitive to ischemia (Mascalchi et al., 2005). The tech- 
nique is based on the use of a pair of counter-balanced, pulsed 
imaging gradients applied in the same direction. The first gradient 
pulse is used to encode protons as a function of position, whereas 
the second is used to remove this encoding. If protons are com- 
pletely static, the counter-balanced gradient pulses have no effect 
and the resultant MRI signal is proportional to the net magneti- 
zation within a voxel. However, if the protons move between the 
application of the gradient pulses, as is the case for water mole- 
cules diffusing in biological tissue, then the spatial encoding is not 
fully removed and is measured as a signal attenuation providing 
diffusion-weighted contrast. If the degree of diffusion weighting is 
manipulated in a series of DWI experiments, simple mathematical 
modeling enables calculation of an apparent diffusion coefficient 
(ADC) image, which has information distinct from that provided 
by proton density or relaxation time parameters. 



anatomy. The inset shows medullary draining veins (blue arrows). (D-F) 
Patient 1 provides structural images, T1 (D), FLAIR (E), andT2-weighted (F), 
with the latter two images showing subcortical white matter disease (yellow 
arrow). (G) a normal perfusion-weighted image obtained by 
contrast-enhanced dynamic susceptibility contrast. 



Within minutes of vessel occlusion, failure of the sodium- 
potassium pump and associated cytotoxic edema leads to an influx 
of water from the extracellular space to the intracellular space. This 
impedes the diffusion of water molecules, as cell membranes act 
as more of a barrier. The restricted diffusion of water molecules 
results in increased signal intensity on DWI. The use of DWI is 
approximately four to five times more sensitive for detecting acute 
stroke than is non-contrast CT (Chalela et al, 2007). 

The combined use of FLAIR and DWI helps to distinguish acute 
from sub-acute and chronic stroke lesions. The FLAIR technique 
(De Coene et al., 1992) is a variant of the historic "inversion recov- 
ery" NMR experiment, in which an RF pulse is used to invert net 
magnetization. During the ensuing Tl recovery, the inverted lon- 
gitudinal magnetization initially shrinks in magnitude, vanishes, 
and then grows toward its equilibrium value once again. For FLAIR 
imaging, and timing parameters are adjusted so that the MRI sig- 
nal is acquired when cerebrospinal fluid is Tl -nulled. Given that 
CSF tends to pool within the infarct zone as time progresses, use 
of FLAIR combined with DWI can improve the identification of 
new lesions near sites of prior ischemic injury, potentially provid- 
ing insight into stroke physiology and subtype. At present there 
is some debate on whether the FLAIR on its own can be used to 
predict the time from symptom onset. Typically, FLAIR imaging 
can detect the presence of ischemia approximately 3 h after stroke 
(Thomalla et al., 2009) . Recent work suggests that using FLAIR rel- 
ative signal intensity ratios from ischemic to contralateral regions 
of interest maybe better than visual inspection alone for predicting 
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Imaging acute and sub-acute stroke 




FIGURE 2 | Patient 2 is an acute stroke with mild impairment 
(NIHSS= 1) yet the diffusion-weighted image (DWI) (A) shows 
restricted diffusion in the left frontal lobe (red arrow) and a large 
vessel obstruction on the x (TOF) angiography (B). (CD) Arterial spin 
labeling images show reduced cerebral blood flow (CBF) and prolonged 
arterial transit time (ATT) (shaded in dark red) (Based on Macintosh et al., 
2010). Patient 3 shows the right parietal lobe stroke lesion onT1 (E) and 
T2-weighted (F) images. 



the time of symptom onset. This approach has good sensitivity but 
poor specificity (Cheng et al, 2013). Time of symptom onset is of 
critical importance given that the only clinically approved ther- 
apy for acute stroke at present, recombinant tissue plasminogen 
activator (rt-PA), requires administration within several hours of 
stroke onset (Taussky et al, 201 1). Therefore, maximizing the use 
of rt-PA places an increasing need on the availability and capability 
of MRI for assessing acute stroke. 

Magnetic resonance angiography provides a non-invasive 
method to screen for vessel occlusions, stenosis, or malformations. 
There are two MRI techniques available for this purpose, with 
each performing slightly less sensitively than the gold standard, 
digital subtraction angiography. Three-dimensional time-of-fiight 
angiography (3DTOF) (Davis et al., 1993) provides a non-invasive 
method of assessing the intracranial circulation. The technique is 
based on the fact that if RF pulses are applied rapidly in succession, 
there is little time for Tl recovery of longitudinal magnetiza- 
tion. Over time, a steady-state of longitudinal magnetization is 
obtained that is lower than the equilibrium magnetization. How- 
ever, if blood flows into the imaging plane it still retains its 
equilibrium value and thus exhibits positive image contrast with 
respect to the background static tissues. This effect enables imag- 
ing of the vascular lumen, but only for through-plane flow. An 
alternative to 3DTOF is contrast-enhanced magnetic resonance 
angiography (CE-MRA) (Runge et al., 1993), typically involving 
the intravenous injection of a paramagnetic contrast agent such 
as Gadolinium Diethyl-Triamine-Penta-Acetic acid (Gd-DTPA). 
This agent enhances Tl relaxation, increasing the signal intensity 
of blood on appropriately Tl -weighted images, with respect to 
the signal from static tissues. This technique provides more con- 
trast than 3DTOF and is much less sensitive to flow dynamics. 



Consequently, CE-MRA can be used more effectively to image 
extracranial and intracranial vessels. If scan time permits, the two 
angiographic methods are complementary and use of both 3DTOF 
and CE-MRA can improve diagnostic ability (Sohn et al., 2003). 
Evidence is also accumulating that CE-MRA can be used to scruti- 
nize "beyond the lumen," improving characterization of vulnerable 
plaque in the carotid arteries and potentially providing a marker 
of stroke risk (Wasserman, 2010). A disadvantage of CE-MRA, 
however, is that Gadolinium contrast agent administration is con- 
traindicated in patients with poor renal function (Morita et al, 
2011). 

Perfusion-weighted imaging is similar to CE-MRA in that a 
bolus of Gadolinium contrast agent is administered intravenously, 
although in this case fast imaging (i.e., EPI or a variant) is 
employed to generate a time series of images to track reductions 
in T2* -weighted signal intensity of tissues as the contrast agent 
travels through the microvasculature. Semi-quantitative perfusion 
maps are obtainable from this examination that estimate cerebral 
blood volume (CBV, the blood volume per unit of brain) , the mean 
transit-time (MTT, the average time required by the bolus of con- 
trast agent to cross the capillary network), the cerebral blood flow 
(CBF, the volume the blood flowing per brain mass and per unit 
of time, mL/100 g tissue/min). When combined with DWI, both 
imaging data sets provide information about the ischemic core and 
about the ischemic penumbra. It was originally thought that the 
difference between the spatial extent of perfusion deficit and DWI 
hyperintensity reflected the ischemic penumbra. It is now known 
that DWI hyperintensity can resolve in many stroke patients 
(Labeyrie et al., 2012), and that the extent of deficits observed 
by PWI do not necessarily reflect the true extent of the penumbra. 
There is some consensus that penumbra imaging can be used in the 
medical decision making, as in the case of rt-PA therapy (Donnan 
et al., 2009), however, this imaging approach appears to have gen- 
erated mixed results in recent endovascular mechanical embolec- 
tomy clinical trials (Tenser et al., 2011; Kidwell et al., 2013). 

Lastly, T2* -weighted imaging has multiple applications in acute 
stroke. These applications stem from the fact that the T2* of blood 
varies with oxygenation content, based on the magnetic suscep- 
tibility characteristics of hemoglobin. Oxygenated hemoglobin is 
diamagnetic, causing a slight decrease in the applied magnetic field 
that is supported within blood, whereas deoxygenated hemoglo- 
bin is paramagnetic, concentrating the magnetic field within red 
blood cells. Abnormal accumulation of deoxygenated blood thus 
provides hypointensity on T2* -weighted images and is an indica- 
tor of vascular pathology. Thus, T2* -weighted images are capable 
of detecting acute hemorrhage, with equivalent accuracy to CT 
(Chalela et al., 2007). Microbleeds, indicative of multiple types of 
micro-angiopathy are detected on T2* -weighted images (Fazekas 
et al., 1999), but not on CT due to insufficient signal contrast and 
spatial resolution. T2* -weighted images are also capable of depict- 
ing hemorrhagic transformations of ischemic infarcts, and provide 
indications for ruling out primary hematoma, as well as depicting 
thrombosed veins or sinuses. In recent years, T2* -weighted appli- 
cations are being enhanced through use of a technique known 
as susceptibility-weighted imaging (SWI; see Figure 1C), which 
provides improved characterization of susceptibility changes in 
the brain microvasculature for angiography, venography, and 
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detection of atherosclerosis and thrombosis (Barnes and Haacke, 
2009). 

MRI IN STROKE RECOVERY 

Despite the continuing technical advances in MRI of acute stroke, 
these developments represent only the "tip of the iceberg" with 
respect to health care management of all stroke patients. Currently, 
only a small fraction of patients are eligible for acute stroke ther- 
apy. The remaining stroke survivors are often left with profound 
deficits, leading to substantial loss in quality of life, economic 
productivity and substantial caregiver burden (Broderick et al., 
1 989) . In the developed world, the problem is worsening as seniors 
become more prevalent due to the "baby boom" generation and 
other factors. 

Although MRI is often viewed as not widely accessible and 
expensive, its versatility and non-invasiveness makes it well suited 
for studies involving stroke patients to understand stroke pathol- 
ogy and the process of recovery. If useful therapeutic applications 
arise from such studies, then MRI could be used potentially as a 
biomarker in the clinic. Used sparingly and appropriately, the costs 
associated with additional MRI examinations in stroke patients 
will be small compared to the costs of caregiving, and any substan- 
tial reduction in the extent of disability due to stroke introduced 
by such new therapies will be very important. In the remainder of 
this review, representative examples are provided that illustrate the 
breadth of the opportunity, and also the challenges that lie ahead. 

MRI FOR NEUROANATOMICAL CORRELATES OF STROKE 
SYMPTOMS 

Given the spectrum of deficits exhibited by stroke patients, the 
variability in location, and extent of ischemic stroke, and the cur- 
rent imperfect understanding of normal brain function in the 
young and in the elderly, it is not surprising that there are a substan- 
tial number of studies that investigate how stroke lesions in specific 
brain locations correlate with behavioral deficit. The use of stroke 
patients as a "knock-out" model of brain function can help to clar- 
ify the role of a given brain regions, although many strokes are not 
sufficiently punctate for definitive conclusions and findings must 
be corroborated through replication and consensus. For stroke 
patients, improved knowledge of how damaged neuroanatomical 
structures relate to behavioral deficits can potentially be important 
for refining cognitive and physical rehabilitation strategies, in the 
long term. 

A simple example illustrating use of MRI in this context involves 
use of high resolution Tl -weighted images to compare the lesion 
location among chronic ambulatory stroke survivors with sym- 
metric gait compared to patients with asymmetric gait (Alexander 
et al., 2009). Individuals with asymmetric gait were more likely 
to have the putamen included in the lesion volume, a structure 
essential for motor control (see Figure 4). 

Such results are increasingly made possible by new develop- 
ments in image processing that facilitate volumetric analysis of 
anatomical MRI data. One example is voxel based morphometry 
(VBM), an approach that involves spatially registering anatomi- 
cal MRI data from a group to a known "atlas" of brain regions, 
and then comparing image volumes across brains at each voxel 
(Ashburner and Friston, 2000). Another approach is measuring 



the cortical thickness by segmenting cortical gray matter from 
CSF, and cortical gray matter from subcortical white matter, based 
on signal intensity differences in high resolution Tl -weighted 
images (Fischl and Dale, 2000). Tissue boundaries can then be 
represented as convoluted surfaces, with point-to-point thickness 
estimates obtained by projecting orthogonally from one surface 
to the other. A recent study used VBM to help classify stroke 
patients that responded positively to epidural motor cortex stimu- 
lation compared to non-responders (Kansagra and Wong, 2008a). 
One study of cortical thickness found that patients with severe 
cerebrovascular steno-occlusive disease who underwent surgical 
revascularization appeared to have replenished the cortex with 5% 
increase in thickness (Smith et al., 2010). Some of the challenges 
of cortical thickness methods in cerebrovascular patients relate to 
the available signal contrast and SNR at 1.5-3.0 Tesla, which may 
be addressed by ultra-high field MRI (Zwanenburg et al., 2012). 
Furthermore, misclassification of tissue, advanced atrophy, or in 
the presence of multiple lesions are additional challenges that may 
necessitate within-patient analysis approaches that avoid the image 
processing step of group template coregistration (Longstreth et al., 
1998). 

DIFFUSION TENSOR IMAGING 

It has been known since the 1990s that there is a preferred direc- 
tion for water diffusion in certain tissues (i.e., the diffusion is 
anisotropic) (Henkelman et al., 1994). The effect is pronounced 
in white matter tracts, with diffusion preferentially occurring 
along the length of axons, and more restricted in the orthogo- 
nal directions. The myelination of axons, the axonal membrane, 
and microtubules within axons are all thought to contribute to the 
diffusion anisotropy. If a series of DWI scans are conducted with 
diffusion weighting in a variety of different orientations, it is pos- 
sible to estimate the direction and magnitude of diffusion and the 
correlations between diffusion in different directions, properties 
that define the anisotropic diffusion tensor and hence lead to the 
term "diffusion tensor imaging" (DTI) (Nucifora et al, 2007). A 
minimum of 7 DWI experiments with different diffusion orienta- 
tions (including one experiment with no diffusion weighting) are 
required to calculate diffusion tensor parameters, although encod- 
ing for more diffusion directions improves data quality. Thus, DTI 
experiments are often time-consuming (approximately 10 min or 
longer), even when EPI techniques are used for spatial encoding. 
The technique can be sensitive to a variety of different image arti- 
facts, including head motion, as the signal attenuation measured in 
DWI experiments does not differentiate well between movement 
of water molecules on a microscopic or a macroscopic scale. 

Data provided by DTI enable calculation of "rotationally invari- 
ant" signal contrast parameters that are theoretically independent 
of which gradient orientations are chosen for diffusion weighting 
(as long as they are well-spaced in 3D). The mean diffusivity para- 
meter provides contrast similar to that of the ADC, by providing 
a weighted average of diffusion in three orthogonal directions. A 
parameter known as the fractional anisotropy (FA) describes the 
extent that there is a preferred diffusion direction within a partic- 
ular voxel (shown in Figure 3A). The FA value for gray matter is 
thus lower than that for white matter, potentially providing a rich 
source of information relating to brain connectivity. 
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FIGURE 3 | Diffusion tensor imaging can be used to visualize the 
fractional anisotropy (FA) (A) in white matter (brighter signal means 
higher FA) or produce the white matter tracts (B). In this case an atlas 
based corticospinal tract is shown in red and overlaid on a coronal 
anatomical image. Functional MRI results shows task-based activation (C) 
and a resting-state network corresponding to the default mode network (D). 



Location of the stroke lesion will strongly influence the poten- 
tial for recovery. For example, from lesion analysis work, it is 
known that a stroke lesion to the internal capsule is particularly 
problematic for motor recovery because the output of motor sig- 
nals converges on this subcortical brain structure (Shelton and 
Reding, 2001). One of the major appeals of DTI for stroke recov- 
ery research is therefore the ability to visualize the integrity of the 
corticospinal tract (shown in Figure 3B). DTI can be used to probe 
the integrity of the corticospinal tract within the first weeks after 
stroke, and early DTI measures are highly correlated with residual 
motor function in the acute and chronic stage (Radlinska et al, 
2010). Thus, DTI measures may help to inform whether a patient 
has the potential to recover a particular behavior or not. 

Diffusion tensor imaging data can also be used for tractogra- 
phy, which is a non-invasive method for mapping the paths of 
white matter fibers in the brain (Ashburner and Friston, 2000). 
Remarkable pictures of white matter tracts have been generated 
in monkeys (Adluru et al., 2012) and in the post mortem human 
brain (Kansagra and Wong, 2008b) although this image quality 
is currently logistically impractical in routine clinical MRI due 
to the long scan times that are required. Part of the reason for 
these long scan times is that the basic diffusion tensor model 
for white matter does not account well for regions of the brain 
that have crossing fibers. More sophisticated models require more 
gradient orientations, and high spatial resolution is paramount. 
This is a promising methodology for use in stroke patients, and 
concerted research efforts are under way to improve MRI hard- 
ware and enhance the clinical practicality of DTI tractography 
(Van Essen et al, 2012). For now, however, FA measures derived 
from DTI remain the primary method of quantifying anisotropic 
diffusion. 

FUNCTIONAL MRI 

Blood oxygenation level-dependent (BOLD) fMRI (Ogawa et al, 
1990; Borresen and Lambert, 2008) has revolutionized neuro- 
science through its capability to measure signal changes associated 
with neuronal activity generated by sensory stimuli, or by behav- 
ioral tasks involving memory, cognition, action, or emotion. Neu- 
rons communicate with glial cells and the nearby microvasculature 
to signal for the delivery of additional blood through vasodilation, 
when neuronal activity levels increase. The neurovascular unit that 
governs such processes represents the ultimate spatial resolution 
achievable with hemodynamic measures of brain activity, in the 
range of hundreds of microns (Menon and Goodyear, 1999). Due 
to contrast to noise ratio considerations, however, fMRI studies are 
typically conducted with voxel dimensions of several millimeters. 
Researchers have determined that the local field potentials arising 
from healthy populations of neurons are strongly correlated with 
BOLD-fMRI signals (Logothetis et al, 2001). 

In healthy individuals with an intact neurovascular unit, the 
BOLD-fMRI response to neuronal activity is characterized by three 
distinct phases: (1) a fast response lasting 1-2 s when there is a 
very small decrease in the BOLD signal, (2) a larger hyperemic 
signal increase caused by the inflow of oxygenated blood, peaking 
approximately 4-5 s after the neural stimulus, and (3) a refractory 
period lasting approximately 10 s during which the BOLD signal 
undershoots and then reaches the baseline level (Norris, 2006). 



The period of hyperemia provides the most robust means to detect 
increases in brain activity and the shape of this hemodynamic 
response is often modeled as a gamma-variate function (Cohen, 
1997). BOLD signal characteristics are a complex function of mul- 
tiple physiological parameters such as CBF, CBV, cerebral rate of 
metabolic oxygen consumption, and hematocrit. The key property 
used for generating fMRI, through which these other parame- 
ter changes are viewed, is the effect of oxygenation status on the 
net magnetic susceptibility of blood (oxyhemoglobin is diamag- 
netic and deoxyhemoglobin is paramagnetic, as described above). 
Thus, the predominant method for measuring BOLD-fMRI signal 
changes at 1.5 Tesla and 3.0 Tesla involves using T2*-weighted EPI 
technique to generate time-series data during periods of neuronal 
activity. The resulting signal changes are typically small, only a few 
percent from signal baseline, partly due to the small volumetric 
fraction of blood in tissues. The upper bound for the BOLD sig- 
nal, which occurs when deoxyhemoglobin concentrations become 
negligible, has been estimated as approximately 6% in gray matter 
at 3T (Bulte et al, 2012). Instances where computed BOLD-fMRI 
signal changes are larger than 6% should be scrutinized rigorously 
for spurious signal contamination, from sources such as larger- 
scale vasculature (Nencka and Rowe, 2007) and head motion. 
Regarding the latter source of artifact, it is unfortunately the case 
that head motions of a few millimeters easily contaminate BOLD- 
fMRI signals. This has long been a problem that affects clinical 
utility of fMRI, although recent technological improvements are 
beginning to have an impact (Rotenberg et al., 2013). 

TASK-BASED BOLD FMRI 

The majority of the fMRI literature includes task-based 
experiments, whereby participants alternate between different 
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FIGURE 4 | (A) Overlaid on the average T1-weighted image is a color map showing the average lesion location as a proportion of stroke patients with normal 
gait. (B) Among stroke patients with asymmetric gait, there was a higher prevalence of stroke lesion that included the putamen and corona radiata. (Based on 

Alexander et al., 2009). 



behavioral or stimulus conditions to induce a measurable BOLD 
signal change. For example, the simplest approach is to alter- 
nate between one specific stimulus/task and a resting condition 
(often visual fixation on a displayed cursor). When the stimu- 
lus/task is very brief (approximately 0.1-4 s), the experimental 
design is described as "event-related," whereas longer duration 
task and control condition periods (typically 15-30 s) are used in 
"block-designs" (shown in Figure 3C). In stroke recovery research, 
the majority of fMRI studies have used a block-design primar- 
ily because the block-design affords greater detection sensitiv- 
ity (Calautti and Baron, 2003). Irrespective of the experimental 
design, the task of interest is undertaken in multiple repetitions to 
improve the ability to detect brain activity under limited SNR con- 
ditions. The majority of fMRI analyses use a simple General Linear 
Model approach to determine which voxels in the brain show sig- 
nal patterns consistent with the task design, while accounting for 
the sluggishness of the BOLD hemodynamic response (Price et al., 
2006). 

The first task-based fMRI experiment involving stroke patients 
dates back to 1997, in which Cramer and colleagues (Cramer et al., 
1997) visualized the brain regions involved when stroke survivors 



moved their affected limb. Since this pioneering work, there have 
been numerous seminal fMRI studies that have helped to charac- 
terize brain activation patterns that indicate positive or negative 
outcomes after stroke. For example, one sign of poor recovery 
is significant activation in the unaffected hemisphere (i.e., con- 
tralateral to the side of the lesion) when performing a unilateral 
task with the affected limb. Another is diffuse activation of brain 
regions during the execution of a motor task, in regions typically 
not observed in healthy individuals performing the same behav- 
ior. These observations are supported by a long list of evidence in 
the fMRI literature, reviewed by Calautti and Baron (2003) from 
cross-sectional and longitudinal fMRI studies. 

One of the challenges of task-based fMRI in stroke recovery 
research is the issue of task performance. It is important to char- 
acterize properly what the patient "does" during fMRI, to rule 
out the possibility that their activation patterns are different from 
normal individuals simply because they did the task differently 
(e.g., speed and extent, applied force, and level of cognitive effort 
when performing a motor task). One strategy to circumvent this 
important confound for motor fMRI is to measure relevant bio- 
physical signals concurrently, such as electromyography (EMG) 
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FIGURE 5 | A task-based fMRI example from a patient with a right side 
stroke affecting their right leg. (A) An activation pattern generated using a 
standard fMRI model and (B) an activation pattern using an alternative 
approach, that is the patient's electrodermal activity time-series data 
measured concurrently. (C) Moving the right affected limb (R:Aff) resulted 
in larger and repeated phasic EDA responses compared to the left 
unaffected (LUnaff) limb, where EDA responses are small and habituate. 
(D)The event-related fMRI paradigm is shown for reference. These results 
suggest moving the affected limb increased arousal and/or effort (Based on 
Macintosh et al., 2008). 



(Macintosh et al, 2007) or electrodermal activity (EDA) (Macin- 
tosh et al, 2008). EMG provides important information regarding 
the timing and intensity of muscle contraction, and potentially the 
muscle groups that are activated, if multichannel recordings are 
undertaken. EMG maybe used to determine if mirror movements 
confound the interpretation of brain activation patterns, when 
unilateral movement is prescribed. Alternatively, EDA recordings 
provide an indirect measure of the autonomic nervous system, 
probing aspects of arousal, emotion, and sense of effort. EDA 
recordings have been used to characterize sense of effort during 
the performance of a motor task with the affected limb post stroke, 
as represented in Figure 5. 

Despite the wealth of information that can be gained in such 
studies, use of task-based BOLD fMRI to study brain activation 
in individuals with cerebrovascular disease is not without con- 
troversy. Recalling that BOLD signal changes are influenced by 
factors such as CBV, CBF, and the cerebral rate of metabolic oxy- 
gen consumption, pathological changes in these parameters that 
affect fMRI signals maybe mistaken for changes in neuronal activ- 
ity. For example, a reduction in baseline CBF, due to chronic 
hypoperfusion will contribute toward a larger BOLD signal if 
neurons remain viable and active in the region. Beyond hemo- 
dynamic parameters, medications, and various other stimulants 
are known to influence the BOLD-fMRI signal (Pattinson et al, 
2009). It is also possible that the temporal characteristics of the 
BOLD response can be affected after stroke, with evidence for this 
that dates back a decade (Carusone et al., 2002). Patients with 
significant extracranial atherosclerotic stenoses or occlusions may 
have a delayed, abnormal BOLD signal response relative to controls 
(Donahue et al., 2010), due to potential factors such as collateral 
vessels or vascular steal. It may be possible to overcome some of 
the limitations of BOLD fMRI in these cases by performing addi- 
tional fMRI experiments that rely on non-BOLD contrast, such 
as using ASL-based methods for CBF-fMRI (see below) or vas- 
cular space occupancy (VASO) CBV-fMRI (Kleim et al, 2002). 
This is an area of functional imaging that needs to be devel- 
oped if we are to improve our understanding of recovery after 
stroke. 

Another challenging issue associated with task-based BOLD- 
fMRI relates to observances of absent BOLD signals despite 
evidence of existing neuronal activity recorded by an alterna- 
tive functional neuroimaging modality. Others have found in 
stroke that the amplitude of the motor-related BOLD response 
is correlated with the capacity of the vessels to dilate, that is 
the cerebrovascular reactivity (Krainik et al., 2005). Rossini et al. 
(2004) used BOLD fMRI, magnetoencephalography (MEG), and 
transcranial doppler (TCD) ultrasound to study chronic stroke 
survivors and found that although MEG-evoked responses to 
somatosensory stimuli were well detected, an equivalent BOLD 
signal was not seen in a subset of patients. In these patients, 
the absence of task-related BOLD signal was suggested to arise 
from impaired autoregulation of cerebral perfusion, secondary 
to cerebrovascular disease. These studies argue for the impor- 
tance of understanding the underlying vascular health of patients 
when conducting task-related BOLD fMRI, and the need to con- 
firm findings through use of multiple functional neuroimaging 
modalities. 



RESTING-STATE BOLD FMRI 

Resting-state fMRI (rs-fMRI) is a new branch of human brain 
mapping that has emerged in recent years (American College of 
Sports Medicine et al, 2010), although the seminal work occurred 
in 1995 (Jezzard and Balaban, 1995). A major advantage of rs-fMRI 
is that it does not require the patient to perform a behavioral task 
as part of the examination. Thus, concerns regarding abnormal 
behavioral performance of patients are removed. Brain activity 
in the resting state is estimated not by task performance, but 
by low frequency synchrony in the baseline fMRI signal between 
different brain regions, usually in the frequency range between 
0.008 and 0.08 Hz. Maps of spatial synchrony can be obtained by 
"seed-voxel" approaches, in which the strength of correlations are 
reported between the signal from one voxel of interest (e.g., in 
primary motor cortex) and the signal from all other voxels; or by 
multivariate approaches such as independent component analysis 
(ICA) (Beckmann et al., 2005). In either the seed voxel or mul- 
tivariate approaches, the resultant large-scale brain networks are 
called resting-state networks (RSNs), which are thought to reflect 
not anatomical connectivity but instead "functional connectivity" 
of the brain (shown in Figure 3D). 

Task-based fMRI has been used for stroke recovery research, 
while rs-fMRI studies have historically have been less common. 
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This trend appears to be changing, particularly as rs-fMRI can 
be conducted for patients with profound impairments, rs-fMRI 
analysis methods have become more widely available (Filippini 
et al., 2009) and test-retest reliability has been established (Ogoh 
and Ainslie, 2009). One new phenomenon that would have impli- 
cations in stroke rs-fMRI research is how reproducible the signals 
are as a function of time, referred to as the temporal stationar- 
ity (Willie and Ainslie, 201 1). In the future, it may be possible to 
find synergies between task and rs-fMRI approaches to advance 
our understanding of activation patterns that explain or predict 
recovery after stroke. 

ARTERIAL SPIN LABELING MRI 

Another MRI technique for assessing brain function and physi- 
ology is ASL. Although discovered and developed contemporane- 
ously, ASL is a complex imaging experiment that historically has 
not been used as widely as BOLD fMRI (Williams et al, 1992). 
Recent technical advances have made ASL more popular. Two 
images are required in ASL; in the first image, RF pulses are used 
to invert magnetization in arterial blood proximal to the imaging 
region of interest. The "tagged" or "labeled" blood water subse- 
quently is allowed to flow into the imaging region for a sufficient 
time to reach the microvasculature, causing a small decrease in the 
measured MRI signal intensity in proportion to the Tl relaxation 
properties of the tagged blood in the microvasculature, and the 
microvasculature volume fraction within a given imaging voxel. 
The second image (the "control") has no effective labeling but 
otherwise near identical MRI signal contrast weighting. The dif- 
ference between the control and tag images provides an image 
that is CBF-weighted. Although there are several parameters that 
must be measured, or more typically assumed, absolute quantifi- 
cation in ASL is relatively straight forward (Wong, 2005; van Osch 
et al, 2009). The ASL signal difference (i.e., control - tag images) 
divided by the initial magnetization (i.e., related to the proton 
density weighting) is directly proportional to perfusion in units of 
mL/100 g/min. 

As for BOLD fMRI, the ASL perfusion signal is limited by the 
small volume fraction of microvascular blood in tissue (1-2% 
of the available signal in gray matter). Low SNR is a limitation 
of ASL and good temporal image stability is required, as many 
images are typically acquired in a time series and then averaged 
together to yield a perfusion image. Refinements in MRI hard- 
ware and acquisition methods have improved ASL image quality 
due to: (1) improved gradient technology for fast imaging; (2) 
availability of higher main magnetic field MR systems (e.g., >3 
Tesla) (Petersen et al., 2010) that lengthen Tl relaxation times and 
enhance labeling efficiency; (3) multi-slice capabilities (Gunther 
et al., 2005) that improve the volume of coverage; and (4) improved 
theoretical models from which quantitative estimates of CBF can 
be derived (Chappell et al., 2010). 

The ASL technique can be designed to map brain activity anal- 
ogous to BOLD signal contrast (Aguirre et al, 2002), or to produce 
baseline CBF images (Alsop and Detre, 1996) as an alternative to 
PWI, avoiding the use of Gadolinium contrast agents. Two possi- 
ble advantages of ASL over BOLD-fMRI are that ( 1 ) the ASL signal 
is closely associated with CBF, whereas the BOLD signal is of much 



more complex signal origins and; (2) the hemodynamic response 
of the ASL signal to neuronal activity has improved spatial resolu- 
tion. The reduced SNR of ASL images in relation to BOLD-fMRI 
images remains a challenge, but the use of ASL to focus on the CBF 
characteristics of normal and abnormal cerebrovasculature is very 
promising. 

The role that ASL can play in imaging acute and chronic stroke 
patients remains open for exploration, with the first clinical ASL 
stroke study dates back over a decade (Chalela et al, 2000). One 
important consideration is the potential use of ASL instead of 
PWI, particularly for patients in which Gd-DTPA is contraindi- 
cated. Thorough comparisons of the two imaging approaches 
are now being conducted in the literature (Siewert et al, 1997; 
Weber et al., 2003; Zaharchuk et al, 2009) and in diagnostic imag- 
ing departments worldwide, as MRI vendors have recently made 
ASL sequences commercially available. PWI and ASL show similar 
results in healthy brain tissue, but some divergence has been noted 
in cerebrovascular diseases. One option that appears clinically fea- 
sible is the inclusion of a 4 to 6-min ASL acquisition to improve the 
quantitative accuracy of the PWI-derived CBF maps (Debette and 
Markus, 2010). Interestingly, among individuals with suspected 
cerebrovascular disease, roughly half of the patients had an abnor- 
mal ASL CBF image while the DSC image was viewed as normal 
(Zaharchuk et al, 2009). 

The ASL method also appears to have an experimental role in 
assessing minor stroke or chronic stroke recovery. Specifically, it is 
important to identify whether there are delays in the time it takes 
the tagged bolus to travel from the labeling plane to the imag- 
ing plane. This ASL transit-time metric has an analogous PWI 
counterpart, and has been referred to as bolus arrival time, the 
arterial arrival time (Macintosh et al., 2010) and, perhaps most 
commonly, the arterial transit-time (ATT or tA) (Wang et al., 
2003). Vascular pathologies and different vascular status can con- 
tribute to a delay or hastening of ASL bolus. If unaccounted for, 
the delay will manifest as underestimated CBF; decreased ATT will 
have the opposite effect. Some clinical ASL studies have incor- 
porated ATT estimates, although longer scan durations are often 
required. Nonetheless, ATT information may be useful in stroke 
diagnosis (Figures 2C,D, patient 2). Future studies should eval- 
uate how ATT can be used to improve understanding of stroke 
outcomes. 

A recent study used ASL to investigate physiological corre- 
lates of dementia at 6 years post-stroke, compared to adults with 
Alzheimer's Disease (AD) and healthy controls (Kuller et al, 2004). 
In this study structural hippocampal volume was the best predictor 
of AD, whereas the ASL was the best predictor of the post-stroke 
adults going on to develop dementia. One methodological aspect 
of this study was the use of a normalization procedure - the global 
gray to white matter CBF ratio - to increase detection sensitivity 
as is relatively common in SPECT and PET studies. The normal- 
ization of ASL data to provide semi-quantitative images may be 
a better alternative to absolute CBF images in clinical scenarios, 
in which vastly different resting CBF levels are observed across 
patients. Others have shown that use of the normalized CBF image 
had higher clinical sensitivity in a study of AD patients when 
compared to using the absolute CBF images (Benar et al., 2002). 
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FUTURE DIRECTIONS FOR MRI RESEARCH IN STROKE RECOVERY 

Given the breadth and versatility of MRI methodology, there are 
numerous interesting research directions that can be undertaken 
to advance the scientific understanding of stroke recovery, ulti- 
mately toward improving health care management of chronic 
stroke patients. In this closing section, three MRI approaches are 
discussed briefly. 

First, similar to the multifaceted acute MRI protocol, the sci- 
entific field would benefit from improved integration of multi- 
parametric MRI protocols for sub-acute and chronic stroke. For 
example, to date there are few examples integrating DTI white 
matter tract data, rs-fMRI functional connectivity data, and task- 
based fMRI data in stroke patients, yet there should be linkages 
between these datasets. It will be important to assess the extent 
that each MRI pulse sequence provides unique, clinically rele- 
vant information; whether all experiments are required in clinical 
research protocols; or whether only the most sensitive indica- 
tors are practical for expediency. An additional area requiring 
expanded multi-parametric MRI involves improving the quan- 
titative interpretation of task-based BOLD-fMRI data. With the 
potential to develop MRI experiments sensitive to parameters 
such as microvascular CBV and hematocrit, it may be possible 
to combine such measures with ASL-derived CBF to model the 
relationship between BOLD-fMRI signal changes and cerebral 
metabolic rate of oxygen consumption (CMR02) (Hoge et al, 
1999). The prospect of non-invasive CMR02 maps would pro- 
vide a good marker of neuronal pathology for stroke patients. 
Again, the ability to perform such multi-parametric assessments 
in a timely fashion remains a challenge. However, this is quite 
important given the heterogeneity of MRI findings across stroke 
patients, and the wealth of BOLD-fMRI data already collected in 
stroke patients that would benefit from more detailed evaluation 
and insight. BOLD-fMRI studies in stroke patients also require 
additional validation through supplemental measurements using 
modalities that do not rely on the neurovascular coupling phe- 
nomenon. For assessing cortical activity, electroencephalography 
is relatively low-cost and widely available. As in the case of fMRI, 
however, electroencephalography datasets are complex and require 
careful interpretation. 

Second, MRI promises to play a role in the interplay between 
genetics and stroke recovery. Recently, a number of fMRI stud- 
ies have been undertaken to assess the effects of various genetic 
polymorphisms on brain activation patterns. As the role of vari- 
ous genetic factors becomes better understood, it is possible that 
genetic profile may provide some insight into the potential for 



recovery after stroke, and that assists selection of the appropriate 
rehabilitation therapy for specific patients (Pearson-Fuhrhop and 
Cramer, 2010). For example, polymorphisms of the apolipopro- 
tein E (APOE) and brain-derived neurotrophic factor (BDNF) 
genes may influence outcome after stroke (Cramer and Procaccio, 
2012). 

A last area that deserves brief mention is the use of MRI to guide 
intervention and rehabilitation. Technology is continually (albeit 
slowly) developing toward expanding the list of options for ther- 
apeutic intervention for stroke patients in the acute phase. One 
example of emerging technology involves use of focused ultra- 
sound. Through use of ultrasound transducer arrays with high 
channel count, it is becoming possible to distribute heat load 
across the skull while adjusting the transmit phase of individ- 
ual elements such that they constructively interfere to generate 
substantial focused ultrasound pressure fields within the brain 
(Clement and Hynynen, 2002). Such systems require detailed mea- 
surements of skull thickness (as provided by CT, for example) 
to adjust transducer phase appropriately, and MRI guidance to 
assess the effects of focused ultrasound treatment. MRI guidance 
is critical not only to detect the effects of ultrasound on anatomi- 
cal structures, but also because various MRI parameters (and the 
Larmor frequency, in particular) are temperature sensitive (Peters 
et al., 1998). Thus, specialized MRI techniques provide a method 
of non-invasive thermometry that can be used to control focused 
ultrasound heat treatments. Completely non-invasive thermal 
ablation treatments are possible, as demonstrated in patients with 
brain tumors (McDannold et al., 2010), although future applica- 
tions to stroke are possible. One such approach involves use of 
MRI guided focused ultrasound not for thermal therapy, but for 
transient disruption of the blood brain barrier that enables tar- 
geted delivery of stem cells for neural repair (Burgess et al., 201 1). 
Another involves "sonothrombolysis," in which the combination of 
MRI guided focused ultrasound and injected microbubble agents 
are potentially used to enhance the ability of tPA to recanalize 
vessels (Ricci et al., 2012). 

To conclude, MRI is a highly versatile imaging modality that 
can be used in numerous different applications to probe the 
anatomy and physiology of the healthy and injured brain. Sub- 
stantial capacity remains for new technical innovation to expand 
the wealth of biophysical information that MRI provides. In the 
coming decades, it is very likely that MRI techniques will play a 
critical role in the development and subsequent clinical imple- 
mentation of new therapeutic strategies to improve outcome after 
stroke. 
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